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We formulate an extended linear a model of a quarkonia nonet and a tetraquark nonet as well as a 
f~*) ■ complex iso-singlet (glueball) by virtue of chiral symmetry SUl{3) x SUr{3) and t/j4(l) symmetry. 

^SJ ' In the linear realization formalism, we study the mass spectra and components of the low-lying 

scalars and pseudo scalars in this model. The mass matrices for physical staes are obtained and the 
j^ , glueball candidates are examined. We find that the model can accommodate the mass spectra of 

low-lying states quite well. Our fits indicate that the most glueball like scalar should be 2 GeV or 
higher while the glueball pseudoscalar is r;(1756). We also examine the parameter region where the 
\^ , lightest iso-scalar /o(600) can be the glueball and quarkonia dominant but find such a parameter 

r\J ■ region may be confronted with the problem of the unbounded vacuum from below. 



PACS numbers: 



I. INTRODUCTION 



The pseudoscalar, vector and axial- vector as well as tensor mesons of light quarks have been well-understood in the 

^ ' naive quark model in terms of the chiral symmetry. Despite of its success, the naive quark model can not explain the 

f"~^ : scalar meson sector, which have the same quantum numbers as the vacuum. There are about 19 states which are twice 

more than the expected qq nonet as in vector and tensor sectors, while the mass and decay pattern of these low-lying 

C , scalars are different from the expectation of the naive quark model. To understand the nature of these scalars has 

^i ' been the focus of recent studies e.g. see Refs. [il-|j| and references therein. 

C^ Among the low- lying scalar mesons, the lightest scalar /o(600) or a attracts a lot of interests. It is widely believed 

^D ' that /o(600) is like the Higgs boson which plays a crucial role in the spontaneous chiral symmetry breaking. Confir- 

^N I mation of existence of the elusive /o(600) from nn scattering processes settles down a controversy last for more than 

a few decades @, Q- The ttK scattering 'q and analysis from D decay D~^ — >■ K^tt^tt' [7| revealed that n should 

also exist. BES II also found such a n like structure in J/^ decays Q. Combined with the well determined sharp 

resonances, i.e. isoscalar /"(980) and isotriplet a(980) from tttt and Trry as well as KK scattering processes, now it is 

j^ accepted in literature that these low- lying scalar mesons (say less than 1 GeV) can be cast into a chiral nonet. The 

C^ , next important issue is what is the nature of this nonet. 

There are a couple of viewpoints on the nature of this nonet. For example, the tetraquark model [9| can explain the 
mass hierarchy and decay pattern of this nonet quite successfully and is further supported from other experimental 
data, like the photon-photon collision data, which prefer the tetraquark interpretation for the lowest scalar meson 
nonet IC^ (where it is demonstrated that /o(980) should be a tetraquark dominant state with great details). An 
alternative interpretation is that this nonet is bound state of the meson- meson molecule ^lli] . In any way, this nonet 
challenges a self-consistent interpretation in the naive quark model. 

Nonetheless, the agreement on the nature of this nonet has not been achieved yet. For example, recently by studying 
TTTT and 77 scatterings, it is found that this particle could have a sizable fraction of glueball [i2;ll3]- The K- matrix 
analysis [ij] suggested that /o(600)/cr should be a glueball dominant state while /o(980) should be a mixture of 
tetraquark and qq. A recent pole analysis with Fade approximation suggests /o (980) might be more like a molecular 
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state M^ . The nature of this nonet is also a focus of lattice study [T^| . For example, the physical state of a and k 
can have sizable tetraquark components, as demonstrated in a recent lattice simulation [l7[. 

The great success of chiral symmetry in understanding the nature of lightest pseudo-scalars motivates us to the 
extended the linear a models 'l8, IS] to study the nature of these low-lying isoscalar and pseudoscalar (say, scalars 
less than 2 GeV. The historic review on scalars above 1 GeV but below 2 GeV can be found in 1]). Such models 
may shed some lights on the nature of light isoscalars, especially on the issue of mixing among glueball, quarkonia 
and tetraquark. One interesting question which can be addressed by such models is that which of these low-lying 
iso-scalar and pseudoscalar states are more glueball-like. The results from lattice simulations suggest /o(1500) or 
/o(1700) could be a glueball rich iso-scalar while 77(1489) can be a glueball rich pseudoscalar [20|. Since the 0"'""'' and 
^ glueball states can significantly mix with quarkonia and tetraquark states of the same quantum numbers, it is 
necessary to include these states and a glueball state in an extended linear a model. 

In this work, we extend our previous work [18| by including a nonet to accommodate the tetraquark states and 
focusing on the mixing of quarkonia and tetraquark. We attempt to address on the issues which states in the 
pseudoscalar and isoscalar sectors are most glueball-like. Comparing with the systematic work shown in [19| and 
the references therein, we extend the linear a model by including a complex singlet field as a pure glueball field and 
introduce a determinant interaction term [21i | instead of the logarithmic term to solve the C/(l)yi problem as in .22j . 
Our model predicts that the isoscalar glueball should be heavier than 2.0 GeV when the pseudoscalar ry(1726) is the 
best glueball candidates. The lowest isoscalar /o(600) is found to be quarkonia dominant state with a considerable 
tetraquark component. 

The rest of the paper is organized as follows. In section 2, we introduce the extended linear a model and derive the 
vacuum conditions for the condensates of quarkonia and glueball fields. In section 3 we present our detailed numerical 
analysis. In section 4 we close our study with discussions and conclusions. An appendix is provided to show the mass 
matrices of isotriplet, isodoublet and isosinglet states. 

II. THE EXTENDED LINEAR a MODEL 

The extended linear a model can be systematically formulated under the symmetry group S'C/ij(3) x SUj^ii) x U a(X)- 
Three types of chiral fields are included: a 3 x 3 matrix field $ which denotes the quarkonia states, a 3 x 3 matrix 
field $' which denotes the tetraquark state, and a complex field Y which denotes the pure glue-ball states and is a 
chiral singlet. The transformation properties of these fields under the chiral symmetry are defined as follows: 

$ -^ Ul'^U^'r, 

$' ^ Ul'^'U^ , (1) 

where Ul.r are group elements of the SU{3)l x SU{3)r symmetry. While the complex field Y is invariant under this 
SUn{3) X SUl{S) transformation. Under the ^/^(l) transformation, each fields is changed by a global phase factor as 
defined below 

Y^e-^'^^^Y. (2) 

Following the convention of the linear sigma model, we express the quarkonia fields, the tetraquark fields, and the 
glueball fields as: 

$ = Ta(f)a = Tai<Ja + ITTa), 
<P' = TaCJ)',=Taia'a+t7r'a), 

y=-i=(yi+z2/2), (3) 

where matrices Ta — ^ are the generators of C/(3) and A^ are the Gell-Mann matrices with Aq — w I Isxs- Fields, 

aa and tTq, a'^ and tt^, and y and a, denote quarkonia, tetraquark and glueball states in the chiral basis, respectively. 
Up to the mass dimension 0{p*) (it is believed they are the most important operators to determine the nature of 
light scalars of ground states), the Lagrangian of our model can include two parts: the symmetry invariant part Cs 
and the symmetry breaking one Csb'- 

C = Cs + Csb, (4) 



where the symmetry invariant part includes the following terms 

- Airr($t$$t$) _ Ai'rr($t'$'$t'$') _ A2Tr($^$$t'$') _ XyiYY^f 

- [A3ea6ce''^^*/*e'$/'" + h-c] + [kYDet{<^) + h.c], (5) 
while the symmetry breaking part includes the following terms 

CsB = [Tr{B.^) + h.c] + [Tr{B'.^') + h.c] + {D.Y + h.c.) - [A™Tr($$^') + h.c]. (6) 

As evident from our Lagrangian, the choice of these operators are not an exhaustive one and many more terms are 
allowed, as demonstrated in [l9|. In spite of that, these terms can be considered as a leading choice when considering 
the number of quark plus antiquark lines at an effective vertex as argued in [24], where the authors had restricted 
the maximum number of quark plus antiquark lines up to 8. Here we relax such a constraint by including two terms 
related with the tetraquark fields. 

At this point we would like make some comments on the interaction or mixing terms in our Lagrangian. We 
have considered direct mixing terms between quarkonia with tetraquark fields and an interaction term between the 
quarkonia with glueball fields while neglected the direct coupling between the tetraquark field and the glueball. 

The mixing between quarkonia and tetraquark is affected by three terms in this work which are quadratic (propor- 
tional to the coupling Am), cubic (proportional to the coupling A3) and quartic (proportional to the coupling A2) in 
effective fields. The quadratic interaction term in Eq. ^ can be considered as an effective mixed mass term which 
violated the U{1)a symmetry, which is a higher order term if we count in terms of number of internal quark plus 
anti-quark lines. But we do not differentiate between quarkonia and tetraquark effective fields on their underlying 
quark content and treat them on the same footing as effective fields. So we retained these two higher order terms in 
our choice of Lagrangian. Among these three mixing terms only the cubic interaction term has been considered in 
the reference 22]. 

The other term which is different from the work of 22] is the last term (proportional to the coupling k) in Eq. ([5]). 
This choice is motivated from the observation of 't Hooft, who introduced the coupling between a scalar spurion field 
to the determinant of the quarkonia field in order to solve the U{1)a problem. Furthermore, the study of Lattice 
QCD and sum rules reveals that the instanton effect plays an important role in shaping the properties of the glueball 
ground state, it is well-justified to consider this spurion field as an effective glueball field, as argued in flSj . 

It is worth mentioning here that we shall also consider the mixing terms like Tr[$t<I>]y*y and YTr{^^^'). For 
simplicity, we drop those terms in this work but will be studied in future. 

Except the explicit symmetry breaking terms, the chiral symmetry and Ua{Y) symmetry are further broken through 
the formation of condensates. Both quarkonic and tetraquarkonic condensates are responsible for the spontaneous 
breaking of the chiral symmetry from SUl{3) x SUr{3) to S't/y(2) whereas gluonic condensate for the U{1)a symmetry. 
The remnant isospin allows us to represent two condensates for quarkonia and teraquark fields each as: vo, vs and 
v'o, v'g. While the gluonic condensate in our theory is labelled as Vy. 

Expanding fields around these vacuum expectation values we get the expression of tree level potential 
V{vq, vg, v'o, v's, Vy) which should be stable under the variation of condensates, i.e.. 



dV{v,,v',,Vy) 



0, i = 0,8 (7) 



d{vi,v'i,Vy) 
The explicit expressions for each equations can be worked out straightforwardly and are given below: 



^ = oo + — y=(2wo - Vs )Vyk ~ vom^ - [^ + vovg - —=)Xi 



-(wgwo Vs - — -^ + ^— + -^—)X2 - Y ^[^vovo - wsi's )A3 - — A,„ (8) 



g^^-bs- ^.8(.o + ^)vyk - vsm^ - {vo -s-^^ + —, 



i— ^-/|- + -^— + ^ eJT' Y 3 (^8"" + ^0V8 + V2W8W8 )A3 - —Ki (9) 



Isospin 


1 = 1 


I-h 


7 = 


PseudoScalars(P=-l) 


{tt,^'} 


{K,K'},{K\K*'} 


{vi,V2,m,V4,V5 } 


Scalars(P=l) 


{a, a }, 


{k, k'}, {k*, k*'}, 


{/l,/2,/3,/4,/5} 



TABLE I: The categorization of scalar and pseudo-scalar states in term of isospin quantum number are demonstrated. States 
in the same category can mix with each other. 
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vo'vs'^ 



Worst's 



2\/2 



)A2 






(10) 



dV 



b,' 



/ 2 , ,^8^t^o' 

1^8 ^71*' + 



t^o^^s' 



^^s^^s' 



wowsfo 



wo^'s^'s 



6V2 



3\/2 



)A2 



— =(V2woW8 + t'8 jAs - (wo y=- + ^)^'8 Al - yA,; 



(11) 



|^ = V2Z. + J^(^ 



^^0^8 



3^/2 



)fc — Vyiriy^ — Vy'^Xy 



(12) 



These five constraints are nonlinear in terms of condensates vq, vs, and Ug, fg, Vy, but are linear in terms of 
couplings. To avoid solving nonlinear equations, in our numerical analysis, we can choose a set of vq, vs-, and v'q^ v'g, 
Vy as input to solve couplings. 

More precisely, in order to guarantee that our values of {wp, Wg} are physically meaningful, we replace them by two 
positive quantities, i.e. {w^,u(,}. The relation between {v'qjV'^} and {«,,«(,} can be found from ^], and are provided 
below as 



V2 1 



1 , V2 , 



(13) 



The quarkonia condenstates vq and v^ are solved out from the decay constants of Pion and Kaon, which are given 
below 



/o 1 /^ 1 

J- i^ 1 ^ a 1^ I 1 '^ • fl 

= (-/I^O - --;=V8)C0S6'^ - ( — Wq - — W8)sm6l^ 



(14) 
(15) 



These two relations between the decay constants and our model parameters can be found by constructing the Noether 
current and utilizing the PACA relations, as demonstrated in [24| . 



III. NUMERICAL ANALYSIS 



Due to the unbroken S'[/(2)v' isospin symmetry, physical scalar and pseudo-scalar states can be categorized into 
three groups with isospin quantum numbers as / = 1 (triplet), ^ (doublet) and 0, respectively. Only bare quarkonia, 
tetraquark and glueball fields with the same isospin quantum number can mix with each other to form physical states. 
Moreover, there is no mixing between scalar and pseudoscalar fields. Thus the chiral singlet glueball field can only 
mix with the isospin singlets of quarkonia and tetraquark fields. 

Using these facts, the physics states below 2 GeV can be tabulated as given in TablelH where the isodoublet {K^ K'} 
is connected with the isodoublet {A'*, A'*'} by charge conjugation. And a similar relation holds for {k, k'} and {k*, 

K*'}. 

For both isotriplet and isodoublet sector, a 2 x 2 mixing matrix can be used to describe the mixing among quarkonia 
and tetraquark states. While for the isosinglet scalar and pesudo-scalar sectors, a 5 x 5 mixing matrix must be used 
to describe such a mixing. To extract the relevant mass matrices, we use the following substitutions ctq — >■ wq + erg. 



Fields 


■K 


n' 


u 


K 


K' 


fK 


Mass (GeV) 


0.14 


1.20-1.40 


0.15 


0.49 


1.46 


0.13 



TABLE 11: The experimental masses and decay constants of a triplet and a doublet are provided. 



Fields 


a 


a' 


K 


Ac' 


Vi 


V2 


V3 


V4 


^75 


/? 


/2° 


/.° 


/I 


/^ 


Mass(GeV) 


0.98 


1.47 


0.80 


1.43 


0.55 


0.96 


1.30 


1.48 


1.76 


0.4-1.2 


0.98 


1.2-1.5 


1.505 


1.72 



TABLE IIL The experimental mass spectra for triplets, doublets and iso-scalars are used to determine the best fit. 

as — >■ Vs + (Ts, (j'q — > Wq -|- (Tg, CTg — >■ v's + o's, and yi — ?> ?^y + yi while assuming that other fields have no vaeuum 
expectation value. These mass matrices are provided in the Appendix. 

It is useful to count the total number of free parameters in our model. There are 15 free parameters as shown in 
our extended linear a model given in Eq. (|4]) . Most of these parameters can be fixed by the input from the isotriplet 
and isodoublet sectors. The only unfixed parameters are related to glueball sector. Meanwhile, five vacuum stability 
conditions further reduce the number of free parameters. Therefore our model is over-constrained by experimental 
data except the glueball sector. Below we describe how to fix our model parameters. 

• The tetraquark vacuum condensates {v' ,w^} are treated as input and are assumed to be positive but smaller 
than 2 GeV. 

• To fix the parameters in the triplet and doublet sector, we use the physical masses and decay constants of {tt, tt'} 
and {K, K'} mesons, as shown in Table [III The mixing angles for isotriplets and isodoublets, which are labelled 
as 9k and 9-^, are treated as input and are restricted to vary in the range {— f, f} (since it is widely believed 
that Pion and Kaon are quarkonia states. We will address this issue in our discussions). Accordingly, four free 
parameters in our model are fixed with another two free parameters are traded off by two mixing angles. We also 
impose the constraints on the trace of mass matrix of isotriplets a and a', i.e. TrlM"^] = J^i^'^a)^^^ ^^'^ the trace 
of mass matrix of isodoublet k and k', i.e. TrlAl"^] = J^i^'^^)^^^ ■ With these constraints, we choose parameters 
{wq, vs, rri^, (jn'^)^, Ai, A2, A3, X'l, Am, k Vy} as solved out from input of isotriplet and isodoublet sectors. 

• In order to further constrain the parameters related to glueball sector, following the method in the reference 
[25| . we consider two broad conditions from the isoscalar pseudo-scalar sector 



Tr[M^\iodei = Tr[M^%,p , 



Det[Mrj'^ 



Model 



Det[M^^]Exp 



(16) 
(17) 



where, M^ is the mass matrix for the isoscalar pseudo-scalar states. So two more free parameters are fixed and 
we choose {fc,77iy -I- XyVy} as solved from these two constraints. Combined with the solution of kvy given in 
the previous step, the parameter vy is solved out. 

• The five vacuum stability conditions given in Eqs. (|51 ll2p can further help to reduce free parameters in our 
model. In practice, we choose the following five free parameters {bo,bs,bQ,b'g^,D} as solved from these five 
equations. 

After these input, there is one free parameter unfixed which is selected as the glueball mass my. By using it as 
input, we can predict the masses of lowest scalar and pseudoscalar as well as their components. 

Considering that there is a large uncertainty in the determination of tt' mass, we choose to vary this mass within 
the range 1.2 — 1.4GeV. We scan 5 mass values in a 50 MeV step within the above specified range starting with 1.2 
GeV. 

To test our model, here we are going to demonstrate, how well it can reproduce the mass spectra of mesons. Besides 
that we are also interested to see what is the composition of these low-lying scalar. 

There is a huge number of possible solutions in our parameter space. We regard that the best solution for the 
parameter set is the one which closely reproduces the mass spectra of scalars close to the experimental measured 
values. The goodness of the solution is on the basis of smallness of the below defined two quantities: the first one is 
Xi as defined in [23| 



13 M, 



the 



Xi 



E 



M,' 



M/"P 



(18) 



and we also consider the second one which is defined by the least ^ method labelled as xi below: 

2 



13 \M, 



the 



X2 



E' 



M,' 



{SMr^Y 



(19) 



where, M^* '=°('^^p) jg the mass of the each member of scalar or pesudo-scalar family calculated from our model 
(experiment), whereas, 5Mi^^^ is the experimental error for each mass. The sum takes into account 5 pseudo scalar 
masses, 4 scalar masses, the masses of two triplets a and a', and the masses of two doublets n and n' . 

We also take into account the decay width of the lowest scalar /o(600) — ?► tttt as a constraint. This decay width at 
the tree level are computed and those solutions which give the decay width between 0.35 — 0.9GeV are regarded as 
reasonable. 

Below we explain how we choose our best solution. For this we apply two types of minimum xi ^nd X2 analysis in 
two stages. Typically xi h^'S a smaller value while X2 bas a larger value due to the small experimental errors for some 
scalars, like the 771 and 772- Firstly, for a given mass value of m-^' we choose the best fit solution for the mass spectra 
for all scalars. We generate more than 8 million random parameter sets of {9-^, 9k, v'v'g\ to find the best fit solution 
which minimizes the Xi- For each value of rriT^' and each set of {6*^, fl/f, «',«(,}, we treat the bare glueball mass my 
a scanning parameter varying from —9 to 9. After considering the constraint the vacuum must be bounded below, 
i.e. the coupling Ay must be positive, we read out the best fit value for ray- Then we determine the best fit for each 
scanned mjr' from 1.2 GeV to 1.4 GeV with 50 MeV interval. The final best fit is chosen from these 5 best fits with 
the minimum Xi/X2 value. 

The different Xi definitions given in Eqs. (|18m9p yield a similar bestfit, which are presented in Table (jIVp . We 
would like to highlight a few features out from it. 1) It is the negative mass parameter ra\ that triggers the chiral 
symmetry breaking. 2) The sign of vy is correlated with the sign of fc, and the sign of k is determined from the mass 
spectra of pseudoscalar sector. 3) The couplings Ai, A'^, Ay are positive which guarantee the potential is bounded 
from below. 4) The values of Ai, A2, and Ay as well as k are large, which demonstrate the non-perturbation nature 
of the model. 5) The value of Am is found to be negative. 



Parameter 


Value 


Parameter 


Value 


9t: (radian) 


-0.604 


Ai' 


8.248 


9k (radian) 


-0.714 


A2 


76.428 


vo (GeV) 


0.074 


A3 (GeV) 


-0.738 


«8 (GeV) 


-0.115 


Ay 


38.327 


vo' (GeV) 


0.203 


k 


-78.15 


vs' (GeV) 


0.126 


A„ (GeV") 


-1.044 


Vy (GeV) 


-0.109 


bo (GeV^) 


-0.085 


my" (GeV") 


3.0 


&8 (GeV') 


-0.161 


m*" (GeV") 


-0.025 


60' {GeV'-") 


0.166 


m^,'" (GeV") 


0.744 


ba' (GeV^) 


0.18 


Ai 


35.465 


D (GeV^) 


-0.265 



TABLE IV: The values of parameters in our fit are shown where the best value of m^/ is found to be m^i = 1.2 GeV. 

Our best fit result favor the case where the percentage of tetra quark component in tt' meson is about 67.7% and 
in K' meson is about 57.2%. When comparing our result with the previous studies, we find that the tetraquark 
component of K'{1A6) in our result is quite low compared to 95% in [201 and 76% in [23- It would be attributed to 
effects of glueball (note: though glueball does not directly mixed with the other fields in doublet sector, the parameters 
like gluonic condensate and the instanton coupling constant do contribute to the mass matrix in the doublet sector, as 
shown up by kvy) or decay widths of these mesons will put a more stringent constraint on this percentage. For n' the 
percentage of the tetraquark component is in qualitative agreement with that of [25] where they found a tetraquark 
percentage of 85%. 

With the parameter set given in Table (|IVI) . the predicted mass spectra and components for the triplet {a, a'} 
and the doublet {k,k'} are provided in Table IVl It is found that the mass spectra of the triplet are close to their 
experimental values but those of the doublet deviate considerably. The a and k' are more tetraquark-like while a' 
and K are more quarkonia like. 

The predicted mass spectra of pseudoscalars and scalars are shown in Tables (|VHIVIII) . We have a few comments 
in order. 1) The pseudo-scalar mass spectra can fit experimental data better than the scalar mass spectra but the 



tt' Mass (GeV) 


Field 


Our Value (GeV) 


quarkonia (%) 


tetraquark (%) 


Experimental Value (GeV) 


1.2 


a 


1.055 


38.14 


61.8 6 


0.98 


a' 


1.417 


61.86 


38.14 


1.47 


K 


1.13 


62.14 


37.86 


0.80 


Ac' 


1.186 


37.86 


62.14 


1.43 



TABLE V: Mass spectra and components for the triplet and doublet sector based on our fit are demonstrated where the best 
value of m^i is found to be m^i = 1.2 GeV. 

mixing pattern of scalar and pseudo-scalar are quite similar. 2) When the mass rrifp — 1.72 GeV is used to our fit, we 
find a solution with Xy < 0. To guarantee the condition that the potential must be bound from below (i.e. Ay > 0), 
we have to keep mso out from our fit, which explains why in the definition of Xi we only sum over the masses of 
4 scalar. It is found that this condition can predict the lightest glueball scalar should be around 2.0 GeV or so, as 
can be read off from Fig. (HId), while the lightest glueball pseudo scalar should be 771. The mass splitting between 
these two glueball states is controlled by parameters vy and Ay and is found to be around 0.15 GeV. When compared 
with the Lattice QCD prediction for the glueball bare mass reported in [27] where the mass is 1.611 GeV, our result 
my — 1.73 GeV is slightly heavier than this prediction. When my — 1.611 GeV is taken, then the predicted mass 
of the lightest glueball is m^s — 2.29 GeV. 3) The lightest scalar /*'(600) is found to be 0.27 GeV or so and is a 
quarkonia dominant state. 



tt' Mass (GeV) 


J^^ = 0-+ 


Our Value (GeV) 


quarkonia (%) 


tetraquark (%) 


glueball (%) 


Experimental Value (GeV) 


1.2 


V5 


1.858 


0.037 


0.001 


99.962 


1.756 ± 0.009 


Vi 


1.380 


75.803 


24.167 


0.03 


1.476 ± 0.004 


V3 


1.291 


26.700 


73.294 


0.006 


1.294 ± 0.004 


V2 


0.907 


15.852 


84.145 


0.003 


0.95766 ± 0.00024 


Vi 


0.595 


81.607 


18.393 


0.0 


0.547853 ± 0.000024 



TABLE VI: Mass spectra and components for the pseudo-scalar mesons based on our fit are shown where the best value of m^ 
is found to be m^i — 1.2 GeV. 



tt' Mass (GeV) 


J^^ = 0++ 


Our Value (GeV) 


quarkonia (%) 


tetraquark (%) 


glueball (%) 


Experimental Value (GeV) 


1.2 


f 


2.09 


0.01 


0.0 


99.99 


- 


f 


1.487 


77.469 


22.53 


0.001 


1.505 ± 0.006 


r 

h 


1.347 


22.177 


77.82 


0.003 


1.2-1.5 


h' 


1.124 


21.561 


78.439 


0.0 


0.980 ± 0.010 


h" 


0.274 


78.784 


21.211 


0.005 


0.4-1.2 



TABLE VII: Mass spectra and components for the scalar mesons based on our fit are shown where the best value of tti.^/ is 
found to be m^i — 1.2 GeV. 

In Figure [TJ we demonstrate the dependence of Xy and /o masses upon the free parameter m^ with the rest of 
parameters are given in Table pVI) . As shown in Fig. ([T^), when m^ is larger than 3.4 GeV^, the Xy becomes 
negative. Then the potential of our model has to confront with the problem of unbounded vacuum from below. In 
the allowed values of my, the masses of f^ ,i = 1, 2, 3, 4 are almost independent of its value, as demonstrated in Fig. 
([T]d). The upper bound of my is determined from the condition Fjo > 0.35 GeV. 



IV. DISCUSSION AND CONCLUSION 



In this work we develop a consistent model for the scalar mesons below 2 GeV and focus on the mixing effects 
to the mass spectra. In our model we have taken into account of the quarkonia, tetraquark and glueball scalar and 
pseudo-scalar fields. Bare fields with the same quantum numbers are allowed to mixed with each other to form the 
physical mesons. In this way our isospin triplet and doublet mesons are composed of quarkonia and tetraquark states 
and the isosinglet mesons are composed with all the three chiral fields. We have presented our prediction from the 
model for the scalar mass spectra on the basis of two x^ methods and found they yield the similar results. 
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FIG. 1: a) The dependence of Ay upon my is demonstrated. A solid circle marker shows the point Ay — 0, which corresponds 
to TTiy — 3.452 {niY = 1.858) . b) The dependence of mass of /o upon my is demonstrated. A vertical line with my is drawn 
to read out the lowest mass rrifS = 1.86. 

JO 



We also investigated the candidates of the glueball dominant states in our model. What is more encouraging is 
the determined value of the bare glueball mass which is treated as a scanning parameter in our study, which is in 
agreement with the Lattice result 27] quite well. The consequence of the uncertainty in the bare glueball mass is also 
discussed in our work. 

When the constraint for 6,^ and 9k to vary from {— f , f } is changed to {—§,§}, we find solutions with m|, < 
and m% > which can accommodate data quite well but is in contradiction to the general belief that Pion and Kaon 
are quarkonia states. It is also found that when the constraint for Ay > is loosen, we can find solutions that the 
lightest scalar can be glueball dominant. 

In order to develop this model, we have dropped quite a few terms and make some assumption in order to carry 
out the numerical computation. We can extend current study to those cases where the interaction terms between 
glueball and tetraquark field are included. There are more than one choices available to define the interaction between 
different fields, for example, our choice of the instanton induced term is different from that in ref. [22| . It would be 
interesting to show the difference between these two different parameterizations. We included the decay widths of the 
/o(600) — > TTTT to constrain our parameter space, albeit to put a tight constraint on our parameter sets we can consider 
more decay widths of all scalars and pseudo-scalars and even should include tttt and ttK scattering data. In order to 
examine whether our model can accommodate all experimental data, a global fit by treating all free parameters on 
the same footirig in our model is necessary. To extend our model by including the tetraquark field as demonstrated in 
the references [2^, [23] to AdS/QCD framework is also straightforward. Following the previous study |23l. I3ll - l33| . we 
can extend our model to study the role played by tetraquark states in the chiral phase transition at finite temperature 
and finite chemical potential. 
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Appendix: Expression for scalar, pesudo-scalar mass matrix and decay ^vidths 

Different elements of the scalar mass matrix: 






1 
3V2 



' ' I ^ '2 1 \ 

Wo W8 + -Us ) A2 



^(^+v)a3 + -^^ 2 



1 fVoVy ^ VgVy^ , ^ 



V2 



(M,2)33 = m$,2 + (v„'^ + „g'2)^^' + („2 ^ ^2)^^ 



(M,2)44 = m$,2 + L'2 _ y2«o'«8' + ^] A/ + (" 



6 3V2 4 



(M,2)55 = my^ + {3Vy^) A 



Y 



1 



^'^-->'+^<^-sS'^-«Vf->= 






(M.2)i3 = - 



Am + l^ivQVo' + VsVs')X2 + 4W -W0A3 = (Ms )3 



(M,2)i4 
(M,2)i5 



(Ms2)23 



1 

2 

1 r 1 



^(i^gt^o' + ^o^'s' - ^)A2 - 2(^^«8)A3] = (M,2)4i 



Vs 



K^-^)fc ^{M,\ 



-(uswo +W0W8 ^)A2 - 2W-U8A3 



(M,2)32 



^ ^ 2 L " ^3 3V2 3V2 2 ^ ' 



2V2 
71" 



(wo + V2w8)A3l = (Ms^) 



(M,2)25 = i [^(«0«8 + ^)fc] = (Af.^)52 



(Ms2)34 






{A'L 



s M3 



(M,2)35 = (Ms2)53 - (M,2)45 = (^1^2)54 = 



where, M"^ „,,„,, = {M,^)iu M^ 



(M, 



s J22, -M o-o'o-o' 



(M, 



s )33, M aa'aa' 



(M 



s M4: 



M' 



yiyi 



M^ 



(Ms 



M" 



(A/fs')l3, M2,„,3. = (Ms")i4, M'„,y, = (Ms")i5, M2,3,„. = (Ms")23, 



(A.l) 

(A.2) 
(A.3) 

(A.4) 
(A.5) 

(A.6) 
(A.7) 
(A.8) 
(A.9) 
(A.IO) 

(A.ll) 

(A.12) 

(A.13) 
(A.14) 

= (Ms2)55, 
M^ / — 



(Ms')24, M^,,y, = (Ms')25, M^„,,,, = (Ms2)34, M\„,y, = (Ms2)35, M^^,^, = (Ms2)45. 

Different elements of the pseudo-scalar mass matrix: 
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^^-2 



Am -4^/--UoA3 



2\/ gt'sAs 



(M^")41 



31 



1 r 1 / 



15 



2 
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(A.15) 

A.16) 

A.17) 

A.18) 
A.19) 

A.20) 
A.21) 
A.22) 
A.23) 

A.24) 

A.25) 

A.26) 

A.27) 
A.28) 



where, M^„^„ = {M„^)n M\,^, = (M/)22 Af^^„'^o' = (Af^^)33 M\,,^^, = (Af^2)44 M\,y^ = (Af^2)55 
Af\„^« = {Mr,')i2 M\„^„, = (Af^2)i3 Af\„^3, = (Af^2)^^ A//\oj,, = (A//„2)i5 Af^,^^. == (A/^2)23 A/^,^^. = 

(A^^2)24 M^^.y, = {Mr,')25 M^,,,,^, = {Mr,')34 M\„,y, = {Mr,')35 M^^,,y, = {Mr,')i^. 



For the decay constant, we have taken the following standard formula: corresponding to the interaction Lagrangian 
^int = GfoTTpTTp (the subscript 'p' denotes the physical pion fields), the decay constant is given by: 



3s/T 



iM 



Vo 



2 . At the tree level \-iM\ 



{A.29) 



G2 



Where, sj is the symmetry factor, which is in our case is i and kf = 

We calculated this coupling constant from our bare Lagrangian following the procedure presented in [30| The explicit 
expression for the coupling constant is given below (where Rs stands for the rotation mass matrix for scalars) : 
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(A.30) 



The expressions of mass matrices for a — a' and tt — tt' mesons are given below: 
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(Maa'^)ii = m$2 _ AiiXi - B11A2 - C11A3 - DiiVyk 

{Maa'^)22 = ^*'' - ^22Al' - B22A2 
(Maa'')i2 = ^ [A™ - A12A2 - S12A3] 

(M^7r'^)ii = m^^ - -AiiAi - B11A2 + C11A3 + DiiVyk 
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[2 4 
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The expressions of mass matrices for n — k' and K — K' mesons are given below: 
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